Observation of individual molecules trapped 
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Abstract 

For the first time, ordered polar molecules confined in monolayer-deep 
rectangular pits produced on an alkali halide surface by electron irradia- 
tion have been resolved at room temperature by non-contact atomic force 
microscopy. Molecules self-assemble in a specific fashion inside pits of 
width smaller than 15 nm. By contrast no ordered aggregates of molecules 
are observed on flat terraces. Conclusions regarding nucleation and or- 
dering mechanisms are drawn. Trapping in pits as small as 2 nm opens a 
route to address single molecules. 

The adsorption of medium-sized (< 100 atoms) functional organic molecules on 
patterned surfaces has become the subject of intensive studies motivated by the 
prospect of hybrid molecular electronic devices PQ. The operation of such devices 
should be governed by the electronic properties of a single or, at most, of a few of the 
constituent molecules. With this goal in mind, numerous groups have used scanning 
tunneling microscopy (STM) to study a variety of organic molecules adsorbed primar- 
ily on metal surfaces |5| pU |U El El Q- Besides atomic-scale spatial resolution, STM 
provides the means for an "accurate placement of molecules in appropriate position 
and orientation to form a device" PQ. However, as discussed in ref.pQ, many hurdles 
remain on the way towards the development and the integration of useful molecular 
electronic devices. In particular, it is desirable to electrically isolate the device; this 
can be achieved using an insulating substrate. The extension of atomic-scale character- 
ization and manipulation studies to surfaces of insulators is therefore of high interest. 
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It is, however, quite challenging because STM can no longer be applied, unless one 
deals with ultrathin insulating films on metals as substrates [S|. A more fundamental 
motivation for such studies is to understand the delicate balance between intermolecu- 
lar and molecule-substrate interactions which determines the mobility and aggregation 
of molecules, as well as their eventual ordering on ionic substrates. 

Since 1995 non-contact atomic force microscopy (nc-AFM), has proven capable 
of yielding images with atomic resolution on a wide variety of insulating surfaces |J]]. 
This technique relies on a micro-fabricated tip at the end of a cantilever excited at 
its fundamental bending eigenfrequency |10j . Upon approaching the surface, the tip 
first senses an attractive force which decreases the resonance frequency. The negative 
frequency shift, A/, varies rapidly with the minimum tip-distance, especially a few A 
away from the surface. When A/ is used for distance control, contrasts down to the 
atomic scale can be achieved in ultrahigh vacuum. 

Recently, the technique has been used to study molecules adsorbed on various 
surfaces [111 . But so far, only few results have been reported concerning molecules ad- 
sorbed intact on insulators, presumably owing to the high mobility of such molecules 
at room temperature [I"2| . Thus, in order to obtain ordered molecular arrangements 
with a definite orientation on insulating surfaces, it is mandatory to find means to 
lower the mobility of the molecules. In the present work, we report a nc-AFM study 
of chloro [subphthalocyaninato] boron(III) (SubPc) molecules confined on a modified 
KBr(OOl) surface in ultra-high vacuum (UHV) at room temperature. For the first 
time, an ordered arrangement of molecules has been directly observed on a nanostruc- 
tured insulator consisting of intentionally created rectangular pits acting as molecular 
traps. 

Cleaved alkali halide single crystals have been among the first insulating surfaces 
showing atomic periodicity and resolution of atomic-scale defects in nc-AFM exper- 
iments |13|. Clean flat surfaces can be rather easily prepared and imaged in UHV. 
Straight edged pits can be created by exposing the surface to an electron beam for 
a short time |14|. The pits are square or rectangular and one monolayer (ML), i.e. 
0.33 nm, deep on KBr(OOl). First evidence for confinement, albeit without molec- 
ular resolution was found for flat PTCDA molecules 1121 . Here, polar cone-shaped 
molecules are used. The chemical composition and structure of SubPc are shown in 
figs. H]a, b) . The central boron atom is sp 3 coordinated to an apex chlorine atom 
and to three isoindole moieties. The properties and potential applications of SubPc 
and related molecules have recently been reviewed |15| . The choice of this molecule 
has initially been motivated by its large dipole moment (5.4 D) pointing away from 
the protruding electronegative CI atom which might form a strong ionic bond to sub- 
strate cations. Density functional theory (DFT) calculations predict a dipole moment 
close to the measured one and an electrostatically fitted effective charge of ~ — 0.4e 
on CI EJ. 

The dynamic force microscope is a home-built instrument placed in a UHV cham- 
ber with a base pressure of about 5 x 1CP 11 mbar |16|. We used a silicon cantilever 
with a supersharp tip T7f (resonance frequency /o = 165, 513 Hz, nominal spring con- 
stant k = 48 N/m, quality factor Q — 30,800). During the experiments, the tip is 
oscillating with a constant amplitude A ~ 14 nm (±10% calibration uncertainty). 

The KBr single crystal was cleaved in air, transferred to the vacuum chamber, 
heated to 393 K and irradiated with the 1 keV electron beam of a LEED electron source 
to produce the pits |18|. SubPc molecules were sublimed in UHV from a resistively 
heated aluminum oxide crucible while the sample was kept at 353 K and the pressure 
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Figure 1: SubPc molecule: (a) Schematic top view of the chemical structure; 
(b) Side view, with characteristic dimensions and atomic species (CI: green, B: 
pink, N: blue, C: gray, H: white; from ref. 0). 



below 10~ mbar. A quartz-microbalance was used to calibrate the deposition rate. 
One ML of molecules was deposited, with an estimated error of about 20% |19|. 

Figure |5Ja) shows a nc-AFM image of the KBr(OOl) surface at a constant fre- 
quency shift, A/ = —10 Hz, after irradiation, but before deposition of the molecules. 
Monatomic steps and regular pits are clearly visible. Atomic resolution was obtained 
on flat (001) terraces at more negative A/; as expected, only one type of ionic species 
was imaged as maxima |18|. Figures|5Jb, c, d) show images of the sample after deposi- 
tion of 1 ML SubPc. Molecules and substrate ions could not both be resolved with the 
same A/, likely due to different force vs. distance characteristics |13|. The wide-area 
image in fig. p\b) provides an overview of the way the molecules are distributed. They 
decorate steps and edges of wide pits, but only pits with a size smaller than 15 nm 
are filled (rectangular bright areas). No molecules are present on the flat terraces. In 
fig-Hi )) maxima attributed to single molecules can be recognized inside two pits and 
along their edges. The arrow indicates a 2 nm-wide pit trapping only 4 or 5 molecules. 
This is consistent with the nucleation mechanism described hereafter. In fig. Hid), 
the trapped molecules appear self-assembled into regular rows oriented +45° (upper 
pit) or —45° (central pit) with respect to the edges. The distance between adjacent 
rows is 1.4 nm and the distance between two consecutive molecules in a row is 1.0 nm. 
The apparent height of the rows is about 0.6 nm with respect to the surrounding 
KBr(OOl) terrace, or equivalently 0.9 nm with respect to the bottom of the pit. The 
molecules are rather well-ordered in the center of each pits, but the structure appears 
somewhat mismatched along their edges. A blurred bright area attributed to mobile 
molecules and a few protruding molecules are also visible on top of the organized layer. 
When attempting to image at a more negative A/, thereby decreasing the tip-sample 
distance, horizontal stripes appeared, implying that material was dragged by the tip 
while scanning. 

Let us consider the possible arrangement of the molecules inside the pits. Fig- 
ure |2I a) shows that the alignment of SubPc's in rows tilted ±45° with respect to the 
[100] direction of the KBr(OOl) surface can match the underlying lattice if the distance 
between two adjacent rows is 3 x b = 1.40 nm and if two consecutive molecules in a 
row are separated by 2 x b = 0.93 nm, in good agreement with the experimentally 
observed spacings. However, this argument alone cannot explain why molecules are 
only found trapped within the smallest pits and not on the terraces. To address this 
question, we first discuss interactions of a single SubPc molecule with the substrate 
and then those between molecules. Intermolecular bonds are strong and no new co- 
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Figure 2: (a) Non-contact AFM image of the KBr(OOl) surface after electron 
irradiation (frame edge: 150 nm, A/ = — 10 Hz). Monatomic steps and pits, 
one layer deep, are visible. The bar indicates the location of the section shown 
in the inset, (b, c, d)- Non-contact AFM images of SubPc molecules on the 
electron-irradiated KBr(OOl) surface, (b) Frame size: 100 nm, A/ = — 5 Hz. 
The molecules decorate steps and fill pits of width smaller than 15 nm. (c) 
Frame size: 27 nm, A/ = —8 Hz. No molecules are observed on the KBr 
terraces. Molecular self- assembled structures are revealed within two pits. The 
arrow indicates a 2 nm-wide pit trapping a small number of molecules (4 or 5). 
The square denotes the scanning area shown in (d). (d) Frame size: 18 nm, 
A/ = —8 Hz. Molecular resolution reveals that the structures are tilted ±45° 
with respect to the pits edges. Images (c, d) were processed using a technique 
which enhances the molecular contrast. 
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Figure 3: (a) Possible schematic arrangements of SubPc molecules inside a pit 
consisting of 10 x 10 x 1 missing ions. The lattice constant of KBr is a — 0.66 nm, 
thus b = a/V2 = 0.467 nm. The distance between two molecules in the diagonal 
row is 2 x b ~ 1 nm, while the distance between rows is 3 x b ~ 1.4 nm, in 
accordance with the experimental data. The CI apex atom is strongly attracted 
to the electrostatic potential well near the corner site (position I), resulting 
in a tilt of the molecular axis. Position II is energetically unfavorable, (b) 
and (c) Side view along the [110] and [TlO] directions showing the approximate 
orientation of the molecule in position I and II, respectively. 

valent bonds are expected to form. We therefore treat each SubPc molecule as rigid 
and interacting only via electrostatic and Van der Waals plus steric interactions. The 
latter are described by Lennard- Jones 6-12 pair potentials, using UFF parameters and 
combination rules |2()j . 

Steps and pit edges expose ions of alternating sign, thereby preserving overall 
charge neutrality. The resulting electrostatic potential <f> therefore oscillates parallel 
to those edges with a period equal to the lattice constant a = 0.66 nm, as it does 
along [100] and [010] directions above a flat terrace. Laplace's equation implies that 
this oscillating component decays exponentially in empty space. The ions exposed 
at the pit corners (see fig. E3 a )) induce an additional component which varies on 
the scale of the pit width, except close to the corners. These features are readily 
apparent in figs.^Ia-c) which display <j) computed along [TlO], [100] and [001] directions 
passing through the first missing K + ion in the upper left corner (point O in fig. |3J a)). 
The potential was computed for an array of point charges ±e at unrelaxed positions 
of the ions, an approximation which is sufficient for the following semi-quantitative 
discussion. Note in particular the pronounced minimum (actually a 3D saddle point) 
in fig. Ufa) caused by the three Br - ions next to point O. 

The rapid variation of <f> over molecular dimensions requires taking into account 
the charge distribution within a SubPc molecule. In order to describe electrostatic 
interactions, it is sufficient to approximate that distribution by partial atomic charges 
fitted to the electrostatic potential of an isolated molecule in a surrounding shell. 
So-called CHelpG charges )21| were obtained from a density functional computation 
reported in ref.[7j. The appreciable charge on the CI atom (—0.42 e) leads to a strong 
attraction to the upper right corner for a molecule in position I in figs. |3*f a.bh but is, 
however, mitigated by the steric repulsion from the neighboring K + ions. Taking into 



5 




[110] axis (nm) [100] axis (nm) [001] axis (nm) 



Figure 4: Electrostatic potential inside the pit displayed in fig. Of a) along lines 
passing through point O in (a) the [110], (b) the [100] and (c) the [001] direc- 
tions, respectively. The potential was calculated by summing up the Coulomb 
contributions of 100 x 100 x 100 alternating electrical charges ±e at the unre- 
laxed positions of the ions and subtracting the contribution of a 10 x 10 x 1 
squared monolayer slab in the center. 

account appropriate UFF radii, we estimate an electrostatic energy gain of about 1 eV 
for the CI atom. Owing to the decay of (f> and to the alternating signs of the charges 
on the more distant B, C and N atoms, their contribution is an order of magnitude 
smaller. Van der Waals attraction favors the adsorption of a SubPc molecule to the 
bottom of the pit via two of its isoindole "feet" and also its trapping in position I. The 
resulting trapping energy is so large compared to fcsT at room temperature that it 
prevents a SubPc molecule from diffusing away. Therefore corner sites which expose 
K+ ions can act as efficient nucleation centers for molecules like SubPc. In the center 
of the pit, however, the modulation of <j> is only 0.3 V. The corresponding corrugation 
of the potential energy surface is dominated by the electrostatic interaction of the 
protruding CI atom and leads to a barrier of about 0.15 eV, i.e. almost 6ksT. This is 
still insufficient to prevent diffusion of a single molecule at room temperature over the 
time needed to image one (about 12 scanlines, each lasting about 1 s). Moreover, the 
bias introduced by the slowly varying (f> component, apparent in figs.2Ia,b), induces 
a drift of a molecule adsorbed inside the pit towards position I or to the equivalent 
site on the opposite side of the diagonal. This implies that molecules trapped near 
opposite corners should be rotated by 180°, but this prediction is difficult to verify. 
Position II in figs.[3Ia,c) is much less favorable for trapping because charges, or atoms, 
which contribute most to electrostatic or Van der Waals attraction are kept away by 
steric repulsion of the feet from the pit edges. 

The previous estimates are qualitatively consistent with our observations, namely 
that no molecules are imaged on flat terraces, while they appear trapped along the 
edges of steps and of pits smaller than 15 nm. Besides, pits as small as 2 nm appear 
filled; judging from figs. 2] this width is sufficient for the required potential profile 
to develop. The close packing along different diagonals inside the two pits visible in 
fig. Hfd) can be rationalized if those pits expose ions of different sign at their bottom 
right corners. 

In an attempt to understand the observed ordering, we calculated the interaction 
between two SubPc molecules, using the same CHelpG charges and UFF parameters. 
The molecules were aligned as in the row schematically depicted in fig. EJ{a) with 
their CI atoms and two feet constrained to be in Van der Waals contact with a plane. 
The intermolecular binding energy, which is largely dominated by the Van der Waals 
component, is only 0.1 eV at the distance 2b, but increases by another 0.2 eV if the 
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molecules are allowed to reach their equilibrium separation, namely 7.25 A~ 1.56, i.e. 
implying a serious mismatch with respect to the substrate corrugation. Therefore the 
observed stable commensurate order inside pits of intermediate size is presumably the 
consequence of a complex interplay between electrostatic and Van der Waals interac- 
tions between the molecules, the bottom and the edges of the pit, which requires more 
detailed calculations. 

In any case, intermolecular interactions favors the relative orientation of molecules 
in adjacent rows shown in fig. EJa). Note that the image in fig. |2Id) reveals irregu- 
larities, in particular along the pit edges. Once a molecule is trapped at position I 
and additional molecules arrange according to the supposed structure, steric repulsion 
prevents perfect matching along the pit edges (c/. black circle in fig.|3{a)). This mis- 
match must introduce some deviations from perfect ordering even in the absence of 
irregularities in the pit structure itself. Compared to fig. EI a), the larger modulation 
of <j> seen in fig. |l[b) implies a sizable diffusion barrier ~ 0.7 eV along the pit edge. 
This phenomenon, which is a specific property of steps on ionic crystal surfaces |12|. 
might explain why molecules trapped along extended steps, e.g. at the edges of wide 
pits, appear poorly ordered. 

In conclusion, the nc-AFM measurements and theoretical estimates indicate that 
molecules with a protruding atom or small moiety carrying an appreciable charge (ow- 
ing to electronegativity differences between constituent atoms) can be trapped inside 
few nanometer-wide pits created by electron irradiation on a KBr(OOl) surface. Keep- 
ing in mind the goals outlined in the introduction, the next challenge is to study the 
trapping of various molecules, "molecular wires" in particular 3: .22 , equipped with 
judiciously positioned polar groups on different nanostructured ionic substrates. In 
that connection, note that monolayer deep pits have been created on several alkali 
halide (001) surfaces and that their average number, size and separation can be con- 
trolled by adjusting the electron dose and the substrate temperature |14| . A further 
challenge would obviously be to achieve intramolecular resolution. For that purpose, 
beyond the difficulty of the immobilization of the molecules, the tip preparation as 
well as the use of small amplitudes |23| are two key issues that are currently under 
investigation. 
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